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ABSTRACT 

Among radio galaxies containing nuclear dust disks, the bipolar jet axis is generally observed to be 
perpendicular to the disk major axis. The FR I radio source 3C 449 is an outlier to this statistical 
majority, as it possesses a nearly parallel jet/disk orientation on the sky. We examine the 600 pc 
dusty disk in this galaxy with images from the Hubble Space Telescope. We find that a 1.6 ^jm/0.7 /im 
colormap of the disk exhibits a twist in its isocolor contours (isochromes) . We model the colormap 
by integrating galactic starlight through an absorptive disk, and find that the anomalous twist in the 
isochromes can be reproduced in the model with a vertically thin, warped disk. The model predicts 
that the disk is nearly perpendicular to the jet axis within 100 pc of the nucleus. We discuss physical 
mechanisms capable of causing such a warp. We show that precessional models or a torque on the 
disk arising from a possible binary black hole in the AGN causes precession on a timescale that is 
too long to account for the predicted disk morphology. However, we estimate that the pressure in 
the X-ray emitting interstellar medium is large enough to perturb the disk, and argue that jet-driven 
anisotropy in the excited ISM may be the cause of the warp. In this way, the warped disk in 3C 449 
may be a new manifestation of feedback from an active galactic nucleus. 
Subject headings: galaxies: active — galaxies: individual (3C 449) — galaxies: ISM 



1. INTRODUCTION 

Imaging campaigns have established that radio-loud 
elliptical galaxies can host gaseous disks in their central 
few hundred parsecs, thought to be feed ing the massive 
10 9 M^ black holes at their centers (e.g.. IFerrarese et ail 
Il99ll Uaffe eTall ITooll IFerrarese fc Fordl I1999D . Visi- 
ble and radio images of normal and radio-loud ellip- 
tical galaxies can be used to probe the relation be- 
tween the orientation of the jets and the morphol- 
ogy of the gas and dust near the nucleus. Ground 
based studies first discovered a connection, finding that 
dusty disks are often nearly perpendicular to radio jets 
jMollcnho ff. Hummell. fc Bender! 19"9"1 iKotanvi fc Ekersl 
11979(1 . These studies supported the widely held expec- 
tation that there is a link between the angular mo- 
mentum of accreting material and the orientation of 
the jet axis. Subsequent surveys of the observed an- 
gular difference on the sky, ^Sjd, between bipolar jet 
and disk major axes in radio galaxies found a statisti- 
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cally significant peak in the distribution at \l/,m ss 90° 
( Moll enhoff. Hummell. fc Benderll992UKotanvi fc Ekersl 
1979t iMartel et al l 120001 Ivan Dokkum fc Framd 1199.1 
de Koff et all 120001: iVerdoes Kleiin et alJll999(K Taking 
3D jet orie ntation into account and correcting for beam- 
ing effects, iSnarks et all l)2000D also showed that these 
disks are generally orthogonal to jet axes. The correla- 
tion c ould be stronger for FR I type sources than FR II 
types (|de Koff et al J 120001 iFanaroff fc RiievllT974]l . and 
as of yet, no strong trend between the jet axis and the 
galactic isophotal ma jor axis is known to exist in low 
redsh ift objects (e.g.. ISansom et aLlll987t IKinnev et alJ 
2000). 

More recent examinations of 3D jet/disk orienta- 
tions have found additional complexities in the distribu- 
tion. T he statistical work bv lVerdoes Kleiin fc de Zeeuwl 
(2005) found that dust lanes are nearly perpendicu- 
lar to jets, whereas dusty disks with ellipsoidal edges 
have a wide range of intri nsic jet/disk orientation an- 
gles. Schmitt et alJ l)2002t) showed that projection ef- 
fects of radio and gaseous structures on the sky must 
be taken into account in any such study of jet/disk ori- 
entations. This work demonstrated that the statistical 
distribution of ^ jd is consistent with a homogeneous 
three dimensional angular distribution that is modified 
by including a cone of avoidance. The cone corresponds 
to an absence of jet axes lying within ~ 13 — 25° of 
the plane of the disk. This may suggest that jet/disk 
orthogonality is not a trend, but simply a manifesta- 
tion of the cone of avoidance in the underlying distribu- 
tion of relative orientations IjSchmitt et al.l 120021 . The 
notion that disks are intrinsically perpendicular to jets 
l|Verdoes Kleiin fc de Zeeuwll2005|) may also be an arti- 
fact of observational biases and small sample sizes, or 
because only a restricted class of objects exhibit jet/disk 
orthogonality. 
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Fig. 1. — a) HST/ WFPC2 0.7 /im image of 3C 449. Here the highly inclined nuclear disk is clearly visible, b) 4.8 GHz VLA image of 
the radio structure of 3C 449. At a distance of 73 Mpc, one arcminute corresponds to ~ 21 kpc. 



An additional complication in the interpretation 
of any ^jd distribution arises from the sometimes 
disturbed morphology of gaseous structures in the 
heart of man y rad i o galaxies (e.g.. IMartel et alJl2000l 
Ide Koff et 3,1.1120(1(1 l^rTPokkurn fc Franxlll 995ft . Such 
studies suggest that many of these disks are not dy- 
namically stable or relaxed in the galactic poten- 
tial. Detailed observations of gaseous structure in in- 
dividual galaxies have also found evidence for non- 
planar or warped geometrie s (e.g.. Ide Koff et alJEoOOt 
van Dokkum & Franx 1995, Fcrrarcsc fc Fordl 11999 m 
NGC 6251 andlQuillen fc Bowerlll999l in M84). Further- 
more, INoel-Storr et alj 1)2005(1 presents an emission-line 
study of gas kinematics in the heart of radio galaxies, 
finding that many such structures are not consistent with 
a thin Keplerian disk model. 

In a non-spherical quiescent galaxy, dust and gas set- 
tle into a galactic symmetry plane on timescales of or- 
der 10-100 times the rotation perio d at a given radius 
l|Steiman-Cameron fc Durisenl Il988). Previous studies 
have shown that the dust in elliptical galaxies is likely 
external in origin, and settles into the galactic equa- 
torial plane after several dynamical timescales follow- 
ing a merger (e.g.. iTran et al.ll200lUde Koff et alJl2000l 
IGoudfrooii fc de Jond 119951 ILauer et all 1200511 . Gas 
that does not reside in the principle symmetry plane 
of a galaxy must have been injected or perturbed on a 
timescale shorter than this. Since the rotation period 
is short near the nucleus (approximately 2 million years 
at 100 pc in a luminous elliptical galaxy), we would ex- 
pect to find the majority of nuclear gas disks residing 
in a symmetry plane of their host galaxies. As such, we 
would not expect a link between the orientations of ~ 10 2 
pc sized gas disks and radio jets. Instead, we gener- 
ally expect to find settled, planar, and quiescent nuclear 
disks in the symmetry plane of the galaxy. Nevertheless, 
disturbed dust morphologies are common in the heart 
of radio galaxies, and the observed correlation between 
jet and disk axes requires a physical connection between 
these structures over either a small range of orientations 
(e.g., the cone of avoidance model) or for a subclass of 



radio galaxies. 

Rclativistic effects are known to perturb the inner ac- 
cretion region, but the spin of a massive Kerr black 
hole can only strongly affect the orientation of a disk 
within about ten thousand times the gravitational ra- 
dius, r g — GMbh/c 2 ~ 10 14 cm, where M^u is the 
mass of the black hole. Outside a few parsecs, Lense- 
Thirring precession and the associated settling into the 
midplane, called the Bardeen-Petters on effect, requires 
longe r than a Hubble time to operate llLense fc Thirrine 
|1918t iBardeen fc Petter sonl 119751 iCanroni fc Abraham 
2004: lKiimar fc Pringlg [l985[) . The precession rate is 
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where a is the dimensionless spin parameter. At large 
radii, the disk is not expected to align with the equatorial 
plane of the black hole. The inner edge of the disk near 
10r g i s the proposed sit e for jet collimation and acceler- 
ation QRees et alJH982f) . and most studies predict that 
the jet sh ould be aligned with the spin axis of the black 
hole (e.g., IBardeen fc Pettersonlll975l iKumar fc Pringle! 
1985) but not necessarily the angular momentum of the 
disk well outside of r g . A massive spinning black hole 
can be described as an angular momentum reservoir that 
varies in momentum, but only very slowly with the in- 
flux of fuel. The rate of exchange of angular momentum 
between the outer disk and black hole is expected to be 
particularly slow (> 10 9 years) in radio galaxies where 
the black holes are massive (~ 1O 9 M ) and the fueling 
rates are l ow (~ 10~ 4 M^ yr _1 , or well below the Edding- 
ton rate) l)ReesNl978|) . Thus, the angular momentum of 
the outer regions of the disk, hundreds of parsecs away 
from the black hole, is not expected to play a role in 
aligning the jet axis. The apparent statistical bias to- 
ward orthogonal jet/disk orientations therefore presents 
a mystery. We expect the angular momentum of the 
very inner regions of the disk to couple with the spin of 
the black hole and the jet axis, but currently there is no 
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consistent mechanism capable of accounting for jet/disk 
alignment on 10 2 pc scales. 

Noteworthy exceptions to the above generalizations 
have been considered by a number of studies. A merger 
of two black holes could change the spin axis of the ac- 
tive galactic nucleus (A GN ) on a very short tim escale 
lIMerritt fc EkerslEffll ILiull2"004l) . ILubow et all lEffll 
showed that when [3 < h/r, bending waves may cause 
the inner disk to be misaligned with the black hole equa- 
torial plane. (Here, /3 is a dimensionless viscocity param- 
eter, h is t he Gaussian scale height, an d r is the radius of 
the di sk) . iNataraian fc Pringlel l| 19981) ; iScheuer & Feilerl 
(1996) showed that the timescale for a black hole to 
align with a disk should depend on the ratio between 
the black hole and disk angular momentum times the 
Lense-Thi rring precession time, rather t han an accretion 
timescale. iQuillen &: BoweJ l)1999llT997|) showed that the 
pressure in the ambient interstellar medium (ISM) is suf- 
ficiently large to perturb a dusty disk, causing it to pre- 
cess about the jet axis. Structure associated with radio 
jets has been clearly revealed in recent obser vations of the 
X-ray emitting gas in radio galaxies (e.g . . iFabian et al.l 
120031 lHardcastle. Worrall. fc Birkinshawl 11998^ . implv- 
ing that there could be a connection between nuclear 
activity and the energetics of the ISM in radio galaxies 
an d galaxy cluste rs. The radiative instability proposed 
bv iPringla ill 9961) could also cause a di sk to warp (see 
also lMalonev. Beeelman fc Pringld lT996) . 

1.1. The nuclear disk of 3C 449 

It is useful to consider radio galaxies for which the 
jet and disk axes are clearly not orthogonal or are nea r 
the cone of avoidance estimated bv lSchmitT "eTaIll|200l . 
Specimens possessing extreme (nearly parallel) jet/disk 
orientations serve as important test cases for physical 
models capable of coupling the energetics of radio and 
gaseous structures. 

To that end, this paper presents a study of the radio 
galaxy 3C 449. The elliptical galaxy UGC 12064 is the 
low redshift (z — 0.0171) host to this FR I type radio 
source, which was imaged with the Hubble Space Tele- 
scope (HST) WF PC2 camera in the F702W filter cen- 
tered at 0.7 /im (Ide Koff et all feOOO). Figure □ (a) re- 
veals a large, ~ 600 pc dusty gas disk with faint spiral 
features, lying in a plane that is is nearly parallel to the 
jet axis, as seen in Fig.n(b). Projected o n the sky, the 
jet ax is is at a position angle of PA ss 9° l)Feretti et alJ 
1999), running nearly north to sou th, whereas th e ma- 
jor axis of the disk has PA « 160° (|de Koff et al.ll2000|) . 
As such, the jet/disk angular difference ^ jd ~ 29°, es- 
tablishing 3C 449 as a counterex ample to the j et /dis k 
orthogonality trend proposed bv Ide Koff et al.l (2000), 
being the only FR I radio galaxy (out of 9) to possess 
du st features with $,;n < 50°. In the sample discussed 
by iVerdoes Kleiin fc de Zeeuwl l)2005l) . 6 out of 33 ra- 
dio galaxies with measured jet/disk position angles have 
^ jd < 50°. For 3C 449, the angle made by the rota- 
tional axis of the disk and the jet axis is greater than 60°, 
placing 3C 44 9 near the cone of avoidance estimated by 
iSchmitt et all l)2002|) . 3C 449 is at a distance of 73 Mpc, 
such that 1" corresponds to 350 pc for a Hubble constant 
of 70 km s _1 Mpc -1 . The maximum angular resolution 
of the HST data is ~ 0".l, corresponding to ~ 35 pc. 
The mass of the black hole residing at the nucleus of 3C 



449 is estimated at 2.5 x 1O 8 M l|Bettoni et al.l [20031. 
Based on measurements of the central core at optical, ra- 
dio and X-ray wavelengths, the active nucleus of 3C 449 
has a bolometric luminosity ~ 10 42 erg s" 1 or ~ 6 x 10~ 4 
of the Eddington luminosity l|Donato et al.l l2004) . 

The recent HST/NICMOS 3CR snapshot survey 
l|Madrid et all 12005]) allows for a more detailed study 
of the dust features in UGC 12064. Because extinction 
from dust is reduced at 1.6 /im compared to visible wave- 
lengths, the NICMOS image provides a less obscured 
view of the underlying stellar surface brightness profile. 
By combining this near-infrared data with the 0.7 /im 
image, we created a colormap which highlights the ab- 
sorptive features and isocolor contours (isochromes) of 
the disk. We describe these features in §2, and discuss 
an observed twist in the isochromes that is inconsistent 
with the absorptive properties of a planar disk. In §3 
we describe our method of creating model images by in- 
tegrating galactic starlight through an absorptive tiltcd- 
ring disk model. We explore non-planar structures for 
the model disk, finding that a warped geometry provides 
a possible explanation for the peculiarities observed in 
the colormap. Our model is used to constrain the mor- 
phology of the dusty disk, and in §4 we discuss physical 
mechanisms which could account for the warp predicted 
by our model. A summary and discussion follows. 

2. OBSERVATIONS 

3C 449 has been observed as part of a near-infrared 
snapshot survey carried out with the NICMOS 2 cam- 
era on boa rd HST, using the F160W broad band filter 
at 1.6 /im l|Madrid et alJl2005|) . After ali gning and scal- 
ing the images, we created a colormap via division of 
this data with the 0.7 nm image taken with WFPC2 and 
previo usly studied bv lMartef etaf] ll200?il) : lde Koff et all 
(2000). PSF differences between the two images were 
ignored, as any variation would only strongly affect re- 
gions of high surface brightness, (e.g. near the nucleus). 
The images and colormap are shown in Figures ^a) , [21 
and ^ respectively. A highly inclined, sharply edged 
disk can be seen in both the 0.7 and 1.6 /im images of 
Fig.^a) and[21and the correspondi ng colorma p (Fig.lU; 
a structure previous ly described bv lMartel et alJ (2000); 
Ide Koff et all <|200f)J) . The 0.7 /mi WFPC2 image reveals 
more filamentary dust structure than the near-infrared 
1.6 /im NICMOS image. We estimate that the ~ 600 
pc disk has an axis ratio of ~ 0.56 near its outer edge, 
with the 3". 45 long major axis lying at position angle 
PA 20°. 

As expected, absorption from the disk is clearly less 
prominent at 1.6 /mi than at 0.7 /im, so the 1.6 /im NIC- 
MOS image was used to fit a Sersic surface brightness 
profile to the galaxy. We use this profile to represent the 
underlying stellar distribution in our models (see §3). 
The isophotes at 1.6 /im in the central region have ellip- 
ticity ranging from 0.1 at a radius of 2" to 0.2 at 13". As 
can be seen from the contours in Fig. El there is no signif- 
icant twist in the isophote orientation, as each isophotc 
is aligned with major axis at PA w —10°. 

T he major axes of nuclear disks in normal ellipti- 
cal llTra n et al.l 120011) and 3C radio elliptical galaxies 
l|Martel et alJl2000D tend to be aligned within 15° of the 
galaxy isophotal major axis, and 3C 449 is no exception. 
The galactic isophotal axis, lying along ~ —10° at radii 
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Fig. 2. — HST/NICMOS 1.6 fJ.vn image of 3C 449. The near side of a dusty torus is seen in absorption on the western side of the nucleus. 
The jet axis is nearly in the plane of the sky at PA = 9°. 



just outside the disk (of radius ~ 2"), is within 10° of 
the major axis of the outer edge of the disk, at ~ —20°. 
This suggests that the outer regions of this disk could be 
dynamically relaxed in the galactic potential. 

A closer look at the colormap in Fig. 0] reveals the 
surprising isochromal properties of the disk. Note how 
it reddens dramatically from north of the nucleus along 
the major axis (at PA ~ 160°) to the southwestern side, 
below the nucleus. Even more surprising is the "integral- 
sign" like twist in the isochromes bordering the south- 
western edge of the nucleus. This is unexpected, as a 
planar absorbing disk should possess elliptical isocolor 
contours that are aligned with the disk major axis and 
exhibit reflective symmetry about the disk minor axis. 
This is not the case in the colormap of 3C 449. 

Note the linear streaks along the disk in the col- 
ormap, which have been described as spiral arms or 
dusty filaments by previous studies (|de Koff et al.l l2000; 
iMartel et al] 12000'). Near the nucleus, linear dust fila- 
ments are tilted at an angle more nearly east-west than 
in the outer disk. At a radius of 0".2 the linear dust fea- 



tures are at a shallow position angle of ~ 110°, whereas 
near the edge of the disk, the features are nearly aligned 
with the disk edge at PA ~ 160°. For a planar disk, one 
would expect tightly wound spiral features to cause dark 
absorption lanes that are approximately parallel to the 
outer disk edge at all radii. However, if the spiral spi- 
ral features are more open (less tightly wound) near the 
nucleus they may cause a shift in the observed position 
angles of the features. 

The isophotes in the 1.6 (im image of Fig.[3]are nearly 
round in front of the disk (to the northeast), so there is 
no evidence that the stellar distribution is lop-sided or 
barred. The disk outer edge has nearly the same major 
axis as the galaxy isophotes. Outside regions where there 
is strong extinction from the disk, we see no significant 
twist in the position angle of the galaxy isophotes as a 
function of radius. Thus, we have no reason to suspect 
that the galaxy is strongly triaxial in its central kpc-scale 
regions A stellar bar with a major axis of 3"-4" would be 
capable of causing a twist in colormap isochromes, but 
a bar that is sufficiently triaxial to account for such a 
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Fig. 3. — HST/NICMOS 1.6 lira image of 3C 449 with contours marking isophotes. The contours highlight every 0.4 magnitude difference 
in brightness, moving outward from the core. The innermost contour marks a brightness of 13.7 /f-band magnitudes per square arcsecond. 
Note that the major axis of each contour is generally aligned with all others, suggesting that the galaxy is unlikely to be triaxial. 



twist would likely need to rotate to be stable. A rotat- 
ing bar generally truncates sharply near its corotation 
radius. The surface brightness profile on the northeast 
side is smooth (and well fit by a Sersic profile), so it is 
unlikely that the stellar distribution is strongly triaxial 
(rotating or otherwise) in the central few arcseconds from 
the nucleus. As such, the twist observed in the colormap 
is most likely caused by the disk itself, rather than a 
projection effect from a strongly triaxial nuclear cusp. 

This motivates the creation of an absorptive numerical 
model for the disk and underlying starlight distribution. 
In §3 we explore warped geometries for the disk in an 
attempt to reproduce the isochrome twist seen in Fig.Q] 
and constrain the geometry of the disk accordingly. 

3. A NUMERICAL WARPED DISK MODEL FOR 3C 449 

We now describe our procedure for creating an absorp- 
tive tilted- ring model of the gas disk of 3C 449. The 
starlight from behind the disk is integrated, multiplied 
by an absorption factor due to the disk and added to the 
integrated starlight from in front of the disk. To create 
such a model, we require a light density function and a 
geometric model for the disk which specifies its opacity 



as a function of position on the sky and its position along 
the line of sight in the galaxy. 
We describe the stellar luminosity density p(s) where 

s = yj (1 — Cg)r 2 + z 2 . The ellipticity of the galaxy, e g , 

was not allowed to vary with radius. Our function for 
p(s) is consistent with the Sersic law, fit to the surface 
brightness pro file of the NICMO S image using the rou- 
tine GALFIT l|Peng et alJ 12002). The galactic surface 
brightness profile is 



I(s) = exp 



-K I - 



(2) 



where the scale length sq and power n are listed in 
Table ^ (for n — 4 a de Vaucouleur law is recov- 
ered). Here b n = 2n — \ + 0,0 °^ 876 , as approximated 
bv lPrueniel fc Simier] l)1997|) . 

The light density from stars (in 3 dimensions) that 
is consis tent with the Sersic surfac e brightness profile is 
given bv lPrugniel fc Simienl l)1997j) . 
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Fig. 4. — Colormap of 3C 449 generated via division of the 1.6 /im HST NICMOS image by the WFPC2 image at 0.7 /im. Note the 
reddening of the disk south of the nucleus. There is approximately a 0.5 magnitude difference between the faintest and brightest portions 
of the colormap. The radio jet axis is at a position angle of 9°, measured north through east. Note the misalignment of the isochromes, 
particularly the "integral-sign" like twist in those bordering the southwestern side of the nucleus. A planar disk should exhibit isochromes 
that are aligned with the disk major axis at all radii, which is clearly not the case for this object. 



where a r is a free parameter well approximated with the 
expansion a r = 1.0 - i^-p + |^f. We checked that this 
density function provides a good match to the surface 
brightness profile by comparing the model light, gener- 
ated from Eq. |3 and integrated along the line of sight, 
with the surface brightness profile at 1.6 fim and the 
Sersic function that was fit to it. 

The dusty disk is assumed to be comprised of a series 
of rings of radius r, each with an orientation described by 
two parameters: a tilt with respect to the line of sight, 
o;(r), and a(r), the position angle from north on the sky. 
We have defined these angles with respect to the line of 
sight, not with respect to the galactic midplane or the 
jet axis. We assume that u)(r) and a(r) are powers of r, 
such that w(r) cx r a " and a(r) cx r a °. The opacity of 
the disk (if viewed face-on) is described by r(r) cx r" T . 
We specify the angles and opacity at two radii, r m i n and 
Tmax ■ We assume that interior and exterior to these radii 
the disk opacity is zero. 

We create two arrays, with coordinates x, y corre- 
sponding to positions on the sky. The first array con- 



tains the opacity of the disk at this position; the second 
contains the location of the disk along the line of sight. 
Our numerical procedure begins with a random selec- 
tion of points r, 9 in the disk. Using the ring orientation 
functions a(r) and uj(r), the position of this point is pro- 
jected to determine x, y on the sky and z along the line 
of sight. The disk opacity is computed from the radius, 
and stored along with the vertical position of the disk in 
the two arrays described above. 

At each position on the sky, light from the galaxy is 
integrated using the stellar luminosity density p(s) from 
— oo to the disk surface (with position stored in the sec- 
ond array). The sum of light from behind the disk is 
multiplied by e~ r at that position (from the other array) 
to take into account absorption by the disk. This flux 
is then added to the integrated star light in front of the 
disk (from +oo to the disk surface). The opacity at 0.7 
/im of the dusty disk is assumed to be 4.16 times that at 
1.6 /x m, consistent with a galactic extinction law (^lathis 
1990). The output of the light density integration pro- 
duces model images at 0.7 /im and 1.6 /im. As was the 
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TABLE 1 

Parameters for Warp Model 



Parameter Value Comments 



Isophotal Parameters 



n 


1.2797 


From GALFIT 


so 


5". 6625 


From GALFIT 


e 9 


0.175 


From GALFIT 


Opacity Parameters 




0".l 


Minimum HST resolution 


rmax 


1".6 


Estimated from Fig. HI 


T~min 


0.04 


Estimated opacity at r m in 


Tmax 


0.09 


Estimated opacity at r m ax 


a T 


1.0 


For r(r) oc r a ^ 


Warp Parameters 


^min 


-10° 


Disk inclination at r m in 


^max 


-23° 


Disk inclination at r max 


Clu 


1.0 


For ui(r) oc r a " 


&min 


-73° 


PA of disk at r m in 


&max 


-14° 


PA of disk at rmax 


a a 


1.0 


For a(r) oc r aa 



Note. — Isoph otal parameters are from GALFIT 
(Pcim ct al. 2002). Angles for a are position angles 
(PA), measured from north through east on the sky 
(counterclockwise), a; is a tilt with respect to the line 
of sight (w = 0° for edge on disk). A negative angle 
for ui means that the southern edge of the disk is the 
near side. The jet axis is assumed to be in the plane 
of the sky and at PA = 9°. 



case for the HST data, we create a model colormap via 
division of these two output images. 

3.1. Comparison of model and galaxy images 

Our goal in creating models is to test the hypothesis 
that the disk of 3C 449 is warped, and to better constrain 
the geometry of the disk. We alter model parameters to 
obtain a qualitative match to the morphology of both 
the colormap and the 1.6 fim image. Our best matching 
model, shown in Fig. |SJ compares the image at 1.6 ixm 
to the model image, and Fig.Elshows the colormap com- 
pared to the model as viewed in 1.6 /im/0.7 /im. The 
parameters for this model are listed in Tabled 

Figures |S] and [B] illustrate that a warped disk model 
is successful in reproducing the appearance of the HST 
images. Specifically, we find that absorption in a non- 
planar disk can give rise to an "integral-sign" like twist 
in isochromes, as is observed in the colormap of 3C 449. 
The model also successfully exhibits absorption along the 
western edge of the disk in the 1.6 /im image. In TJ-band, 
light from behind the disk is absorbed by the southwest- 
ern (nearest) side, but is visible past the edge. The same 
is true for the representative model, which differs from 
the NICMOS image in that the far side of the disk is 
also visible and appears larger. This is due to extinc- 
tion effects leading to obscuration of the disk at 1.6 /im. 
The colormap reveals the underlying morphology of the 
absorptive dust features, so it was used as the basis for 
shaping our model disk. Because our model lacks az- 
imuthal (non-radial) density variations or spiral surface 
features, it cannot reproduce the linear streaks seen in 



absorption in the 1.6 fim NICMOS image. 

Fig. El is more convincing, as it shows that our model 
can reproduce the "integral sign" twist in the isocolor 
contours near the nucleus of the colormap, as well as 
match the shape of the absorptive features quite well, 
particularly on the southwestern edge of the disk. 

Certainty and degeneracy in model parameters is of 
key importance in judging the qualitative success of the 
model. We first list the parameters that were held fixed 
and then describe those that we varied to find a model 
that best matched the images. The isophotal parame- 
ters n and sq were fixed since they are consistent with 
the Sersic law fit to the surface brightness profile of the 
galaxy. We fixed the galaxy cllipticity e g at 0.175, as this 
produces model isophotes that are sufficiently oblate to 
be consistent with the isophotal ellipticity at 1.6 /xm. 
We set the maximum radius of the disk to match that 
of the outer disk edge, r max = 1".6. Likewise the tilt, 
u m ax — —23° (w.r.t. the line of sight), and position an- 
gle, a max = —14°, of the outermost ring (at radius r max ) 
were fixed and set to be consistent with the location of 
the disk's outer edge. The inner model ring was set with 
r min = 0".l, corresponding to the maximum angular res- 
olution of HST. We cannot probe the geometry of the 
disk at radii smaller than this, so this radius does not 
describe an inner edge for the disk. Rather, it is used to 
provide boundary conditions for the positional functions 
a(r) and w(r), and for the opacity function r(r). 

The model is strongly sensitive to changes in the op- 
tical depth of the disk, which is set by the parameters 
T m i n and T max at the innermost and outermost rings, 
and with an exponent, a T . Such sensitivity allows us 
to constrain the opacity parameters by incrementally al- 
tering T m i n and T max until the isophotes of the model 
colormap match the images. Too high an opacity leads 
to overly strong absorptive features not observed in the 
NICMOS image, and opacity endpoints that are too low 
have the opposite effect. The best match for our model 
was in the range of i?-band opacity r m .;„ = 0.04 ± 0.02 
to T m ax = 0.09 ± 0.02. This physically reasonable R- 
band opacity range generates isochromes that match 
those seen in the HST colormap. Still, there is inher- 
ent uncertainty in our values for these parameters, as 
their effects on the absorptive features differ depending 
on the shape of the disk. A disk of constant density has 
a higher opacity when viewed at high inclination than 
when viewed face-on. Thus, there is a larger range of 
values over which the inner ring opacity will produce an 
"integral-sign" twist. This, along with the orientation of 
the inner ring, is the greatest source of uncertainty in our 
model. 

The model is not as sensitive to the exponent a T in 
the power law r(r) oc r° r . While opacity and ring ori- 
entation effect the structure of the model isochromes, a T 
effects how steeply an absorptive feature tapers off as a 
function of radius. This dependence is minimal, in that 
a model with a T = 0.5 and a model with a T = 1.5 will 
produce similar results, both of which closely resemble 
the data if the positional parameters and opacity bound- 
ary conditions are set appropriately. As such, a T is only 
loosely constrained, and we set it to a T = 1.0 as it gener- 
ates absorptive features in the model colormap that are 
consistent with those seen in the data. 

While the strength of the absorptive features is mostly 
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Fig. 5. — Comparison of a) 1.6 fim NICMOS image with b) the 1.6 (im integrated light warped disk model. Note the similar absorptive 
properties of both the disk and its numerical representation. The disk in the model image appears larger than that in the NICMOS data 
as extinction from the galaxy masks the outer regions of the disk at 1.6 (im. In our model, we do not take into account extinction from 
the surrounding galaxy. The integrated light density function used in creating the model is based upon isophotal analysis of the NICMOS 
image. The i?-band opacity drop-off ~ 0.05 follows a simple power law in an attempt to obtain a qualitative match to the data (see §3). 
Model parameters are listed in Table HI 





Fig. 6. — Comparison of a) colormap of 3C 449, generated via division of the 1.6 ^m and 0.7 fim HST images, with b) the model 
colormap. Note the successful reproduction of the isochrome twist to the north of the nucleus. Our model assumes that the disk is warped 
but featureless, so the spiral features are not seen in the model. Model parameters are listed in Table HI 



dependent upon r(r), isochrome shape most strongly de- 
pends on the structure of the disk. Qualitatively, one 
can achieve similar outer disk appearances by raising a 
and decreasing lu (or vice- versa). Being concerned only 
with the appearance of the disk, the small range over 
which a and w are degenerate is of little relevance. The 
orientation of the innermost ring is of particular inter- 
est, since it is located only ~ 30 pc from the nucleus. 
The best matching model exhibits a dramatic warp, with 
the inner ring at a radius of 0".l at a position angle of 
a m i n = PA = —73°. We estimate the best fit for the 
tilt of the inner ring is at Lo m in = —10° w.r.t. the line of 
sight. 

Remarkably, the best matching model is such that the 
inner regions of the disk are nearly perpendicular to the 
jet axis. Assuming the jet axis is orthogonal to the line 



of sight (e.g., in the plane of the sky) and at PA = 
9°, the axis of symmetry of the innermost model ring 
makes a 3D cone angle with the jet axis of approximately 
10°, meaning the inner regions of the disk are within 
10° of being orthogonal to the jet axis. The uncertainty 
brought about by the codependence of a and ui on disk 
appearance is small enough such that we can be fairly 
certain of this value, as all "good-fit" models have inner 
rings nearly perpendicular to the jet axis. Consequently, 
our model shows that 3C 449 may actually conform to 
the trend proposed by previous studies suggesting that 
dusty disk major axes are preferenti ally oriented nearly 
perpendicular to jet axes on the skv flde Koff et al.l l2f 
Ivan Dokkum fc FranxHl995t iSparks et alJl2000[) . 
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3.2. Comparison of warp model to observations at other 
wavelengths 

The success of the non-planar model in reproducing the 
observed isochrome twist in 3C 449 leads us to suspect 
that the dust disk is warped between a radius of ~ 50 pc 
(the innermost radius we can easily resolve) and 600 pc 
(the disk outer edge). We now discuss the geometry and 
orientation of the dusty disk in context with observations 
o f the nuclear region at other wavelengths. 

iMartel et all l)2000|) presented narrow band images of 
3C 449 in Ha+N[II], showing stronger emission to the 
north of t he nucleus. Likew ise, the near-UV image pre- 
sented bv I Allen et al.l l)2002T) showed similar morphology 
with brighter near-UV emission to the north of the nu- 
cleus. The UV image exhibits arm-like features that are 
bounded by dusty absorptive strea ks also present in the 
visible band image (as discussed bv lAllen et alJ l20021. 

Because extinction is deep along the disk edge to the 
west of the nucleus, we infer that the western side of the 
disk is nearer to us than the eastern side. Studies of the 
radio emission have suggested that the jets are nearly 
in the plane of the sky, at approximately 70° from the 
line of sight, with the sout hern side nearer the observer 
ijGower fc Hutchingslll982]) . As such, we expect that the 
northern jet is in front of the disk and the southern jet 
is behind. The UV emiss ion might be exci ted by star 
formation in the disk (e.g.. lO'Dea et alJl2001l in the case 
of 3C 236), or it might be r elated to shocks or illumina - 
tion along the jet axis (e.g., Solrzano -Iarrea et al.l l'2002'l. 
Since bright UV or Ha+ [Nil] emission is not seen on the 
southern side or along the western disk edge, the emis- 
sion probably does not arise from material embedded in 
the disk. If the emission is from a large diffuse region 
(possibly wide cones) above and below the disk and ori- 
ented approximately along the jet axis, the emission from 
southern side should su ffer more extinction f rom the disk. 
Previous studies (e.g., iQuillen et alJ 11999^ have found 
that structure in an ionization cone can be related to the 
underlying dust distribution in the galaxy. However, the 
UV emission from 3C 449 is brighter outside the dusty 
features than on top of them, though UV emission is 
also brighter in proximity to dust features, so the pres- 
ence of dense material may somehow contribute to the 
excitation. The morphology of the UV emission would 
be difficult to account for with a simple model, but exci- 
tation along a wide cone associated with the jet axis and 
coupled with extinction from the dusty disk might be 
consistent with the observed morphology, so long as the 
UV emission were stronger in proximity to dust features. 
The UV emission also exhibits a deficit (or perhaps a 
shadow) just north of the nucleus. Our model does not 
include material that could absorb UV emission in this 
location, but a warped disk that continues to twist past 
our observational limit at r m i n could possibly account 
for such a shadow. 

We can consider the possibility that blue continuum 
emission associated with the UV structure could have 
reduced the redness of color to the north of the nucleus, 
accounting for the observed isochromal twist in the col- 
ormap. If this were so, we would predict that the size and 
location of features seen in the UV image would corre- 
spond to blue regions in our color map. The colormap is 
quite smooth, wherease the UV emission is clumpy. Fur- 



thermore, bright UV emission is offset from the spiral 
features seen in absorption. This would make it impossi- 
ble for a blue continuum associated with the UV emission 
to cancel out the extinction associated with dust lanes. 

3.3. Comparison to Spectra 

Our warp model predicts that the inner disk (at r ~ 
100 pc or 0".3) is nearly edge-on and oriented almost di- 
rectly east-west, whereas the outer disk at r ~ 600pc is 
oriented with major axis oriented approximately north- 
south. We can compare this prediction to spectroscopic 
observati ons that measure the m ean velocity of the ion- 
ized gas. INoel-Storr et alJ l|2003() obtained spectra with 
STIS on board HST along 3 slits oriented at a position 
angle of 170.6° (approximately north-south along the ma- 
jor axis of the outer disk). Mean velocities measured 
alon g these slits are p resented in Figure 22 and Table 
26 of lNoel-Storr et all (2003). The mean velocity - 0".5 
north of the nucleus is ~ 400 km/s above that ~ 0".5 
south of the nucleus. The observed velocity profile along 
the outer disk major axis is consistent with a rotating 
disk that is approaching the observer on the northern 
side and receding on the southern side. The warp model 
predicts that the eastern side of the inner disk should 
have radial velocity with the same sign as the southern 
side of the outer disk. Because the southern side of the 
outer disk is receding from us, our model implies that the 
eastern side of the inner disk should be receding, while 
the western side should be approaching the observer. 

Mean velocities were also measured in the two slits 
±0."1 from the nucleus and parallel to the central (nu- 
clear) slit. Velocities along the east-west direction can 
be used to test the sign of the rotation predicted by 
the model. The velocities on either side (east and west) 
of the nucleus are lower and higher (respectively) than 
the systemic velocity, suggestin g that the velocity field 
is twisted. East of the nucleus. INoel-Storr et al.l (|2003f) 
found the velocity to be sub-systemic. To the west, a 
velocity consistent with that observed along the major 
axis south of the nucleus was measured. This would be 
consistent with a disk that is rotating such that southern 
and eastern sides of the outer disk are receding from the 
observer. We find that the sign of the velocities along 
the east-west direction is consistent with the prediction 
of our warped disk model. 

4. MECHANISMS FOR CAUSING THE WARP 

We have shown that a warped disk could account for 
the twist in the isochromes seen in the dusty disk of 3C 
449, and that our model is consistent with observations 
at other wavelengths. With an outer edge at r ~ 600 
pc, the disk is approximately in the midplane of the 
galaxy and tilted about 30° away from the jet axis. If our 
model is correct, at smaller radii (~ 100 pc), the disk is 
nearly orthogonal to the jet. IVerdoes Kleiin fc de Zeeuwl 
( 2005) outlined two scenarios that could account for such 
a jet/disk orientation. 

1. The disk is initially planar and sets the jet axis. At 
later times the disk is affected by the torque from 
the galaxy and becomes aligned with the galaxy 
isophotes. 

2. Another force, possibly aligned with the jets, 
pushes on the disk, twisting it at small radii. 
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Might the morphology of the disk in 3C 449 discrimi- 
nate between these two scenarios? We can consider a gas 
disk to be comprised of rings moving in nearly circular 
motion. When the galaxy is non-spherical, it exerts a 
torque on each ring causing precession about a galactic 
axis of symmetry. The precession rate is a ~ f2e$ where 
Q is the angular rotation rate of a particle in a circu- 
lar orbit and e$ is the ellipticity of the galactic gravi- 
tational potential. (Here a refers to a precession angle 
which is measured with respect to a principal plane of 
the galaxy). The precession rate is primarily set by f2, 
and so depends on the rotation period at a given ra- 
dius, P = 2tt/Q. Since the rotation period is shorter 
at smaller radii, the inner disk should be more twisted 
than the outer disk, unless the galaxy ellipticity drops 
rapidly near its inner regions. Such a drop in elliptic- 
ity is not seen in the isophotes of Fig. [31 though the 
presence of the disk itself does not make this absolutely 
certain. As the disk becomes increasingly twisted, the 
inner region can settle into the gal actic midplane (e.g., 
iSteiman-Cameron fe Durisenlll988j) with the outer disk 
warped and misaligned with the galaxy. Instead, we find 
that the outer part of the galaxy is settled into the mid- 
plane and the inner part deviates from or is inclined with 
respect to this plane. This is not consistent with the first 
scenario discussed above. The morphology of the disk in- 
stead suggests that it was once relaxed in the symmetry 
plane and subsequently perturbed by a force from the 
nucleus, pushing the inner disk out of alignment on a 
small, ~ 100 pc scale. 

We now consider candidate scenarios which could ac- 
count for such a warp: 

1. Precession in the inner disk has been induced by a 
binary massive black hole. 

2. Radiation pressure from the central source has 
twisted the disk. The inner disk may have become 
warpe d via a self-induced radiation-driven mecha- 
nism (Pringle 1996). 

3. There is strong coupling between the outer disk 
and the inner disk, aligned by the Bardeen- 
Peter son effect with the spin of the black 
hole ijNataraian fc Pringle] Il998t iScheuer fc Feilerl 
Il99l . 

4. A jet-excited ambient interstellar medium has per- 
turbed the disk. This would be a form of feedback 
from an AGN. 

We discuss these scenarios below. 

4.1. Precession due to a binary black hole 

It is useful to summarize bu l ge an d black hole proper- 
ties for 3C 449. iBettoni et all |2003) have listed a bulge 
velocity dispersion of 225 km and estimated a black 
hole mass of 2.5 x 10 8 Af©. A transition radius (e.g., the 
sphere of influence) where the gravitational force from 
the nuclear black hole dominates that from the bulge is 
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Note that the transition radius can be compared to the 
semi-major axis for a possible black hole binary which 
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where /x is the reduced mass of the black hole binary and 
Mhh now represents the sum of the two black hole masses 
ijCaproni fc Abraham! 12004). For an equal mass binary 
o-hard — »"t/4 ~ a few pc for 3C 449. 

We consider a disk exterior to a binary black hole in 
a circular orbit with semi-major axis db and reduced 
mass \x. At radius r > a, the binary black hole exerts a 
quadrupolar perturbation on the gravitational potential 
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The vertical oscillation frequency is estimated by differ- 
entiating $2 twice with respect to z, and adding it to 
that from the spherical component from the galaxy, 
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where M(r) is the total mass within radius r. The pre- 
cession rate is 
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Relating M(r) to the velocity of a particle in a circular 
orbit, we have M(r) = v 2 r/G. Outside the black hole 
sphere of influence, we expect M(r) oc r, and fl cx r _1 , 
consistent with a flat rotation curve, such that a cx r~ 4 . 
The precession period is then 



10 10 yr 



a b 



lpc 



10 8 M ( 

-2 



100pc 



320km/s 



(9) 



(For an isothermal sphere v c ~ v^cr*.) 

Such a long timescale and steep dependence on radius 
presents a problem for this scenario. For a binary black 
hole semi-major axis larger than a few pc, dynamical 
friction can bring the black holes together on a timescale 
of order a few times the rotation period. Only if the black 
hole binary is hard (a& < a few pc) can the black hole 
separation shrink slowl y, within a timescale on the or der 
of a Hubble time (e.g JMilosavlievic fc Merrittll2003j) . If 
the black hole binary is not hard, then it would spiral 
inward too fast to create the warp. We find that the 
estimated precession timescale is too long to account for 
the warp in the disk. 

One might consider alternative models where the disk 
precessed while a secondary black hole spiraled inward. 
Since the black hole orbit need not be planar, it could 
account for a highly warped geometry in the disk. In 
this case, however, one would expect to see disturbed 
material on both large and small scales. Other than the 
active nucleus, the galaxy appears dynamically relaxed, 
and there is little evidence for a recent galactic merger 
in its isophotal structure (see Fig.[3J). 

4.2. Radiative Warp Instability 
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Here we explore the proposal b y IPringld (Q.996); 
IMalonev. Begelman fc Pringlel l)1996l) that the absorp- 
tion of radiation from a central source can cause a sur- 
rounding disk to precess and warp. The precession rate 
can be estimated from the constant torque solution of 
Maloncv. Bege lman fc Pringlel l)1996f ) (their equation 8) , 
giving a precession period of 
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where E is the disk mass surface density, and h\, i is the 
bolometric luminosity of the central source. Inserting 
values appropriate for 3C 449, we estimate a precession 
period 
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The values used are based on the following estimates: 
The radio core luminosity at 5 GHz is lower than that 
estimated from an unresolved source in the nucleus at 
optical wavelengths, w i th a luminosity of ~ 10 40 8 erg/s 
(JChiaberge et all 119991 ICapetti et aTll2002|) . and at X- 
ray wavelengths wi th a luminosity of ~ 10 41 erg/s 
ijDonato et alJl2004jl . The total bolometric luminosity 
from the core is estimated to be 10 times that in X- 
ray, or L bol ~ 10 42 erg/s llDonato et al.ll2004|) . As- 
suming constant i?-band opacity for a f ace on disk o f 
0.065, we use the extinction coefficients in lMathisI (1990) 
and estimate the surface mass density in the disk to 
be E « 1.3-M© pc~ 2 . The circular velocity v c ~ 
33 0km/s is estimated f rom the velocity dispersion listed 
bv lBettoni et alJ 1)2003!) . assuming an isothermal sphere. 

A quarter precession period would be long enough to 
account for the extent of the warp in 3C 449, and this 
would reduce the above timescale by a factor of a few. 
Even so, this reduced timescale is still very long, so this 
model is not likely responsible for the warp. The preces- 
sion rate might have been faster if the luminosity were 
significantly higher in the past. However, there should 
also be a corresponding increase in the surface density 
of the disk, consistent with the higher rate of accretion. 
We conclude that precession due to radiation pressure 
from a central source is probably too slow to generate 
the warp. The growth time for the radiative warp in- 
stability is approximately the same as the above period 
divided by the order of the mode grown. However, the 
observed warp is large-scale (with wavevector similar to 
the radius), so a higher order mode is unlikely to account 
for it. This finding is cons istent with a similar estimate 
bv lFerrarese fc Fordl ©99), for the warped disk in NGC 
6251. 

4.3. Quick Disk/Black Hole Alignment 

iScheuer fc Feilerl l(T996T) : iNataraian fc Pringlel lfl99l 
showed that a black hole could align with a disk on a 
timescale set by the ratio of the black hole and disk an- 
gular momentum times the Lense-Thirring precession pe- 
riod. The Lense-Thirring precession period is oc r 3 and 
the disk angular momentum is oc r, so the timescale for 
alignment becomes very long at distant radii. The ratio 



of the black hole angular momentum to the disk angular 
momentum is not small enough to reduce the alignment 
timescale at 100 pc to a level that would allow align- 
ment to take place within a Hubble time (see the Lense- 
Thirring precession rate given above in Equation We 
can therefore reject this scenario based on these simple 
timescale estimates. 

4.4. Disk-ISM interaction 

If the pressure gradient from the ambient ISM across 
the dis k is large, it can overcome th e torque from the 
galaxy (Gu illen fc BowedlTiM Il997|) . The pressure re- 
quired to keep a disk of surface mass density E from 
precessing in a galaxy of ellipticity e$ is 



P T >2 x 10~ u dynes cm" 



2 ( e * 



330km/sy \ r 



/100pc 



0.05/ \lM & /pc 2 



cos(6g)sm{0 g ) (12) 



where 9 g is t he disk inclination angle (with respect to the 
g alaxy axis) iQuillen fc Bowerlll999T) 

IHardcastle. Worrall. fc Birkinshawl ©98) observed 
that the X-ray emission in 3C 449 on arcminute scales 
was anisotropic, with isophotes aligned orthogonally to 
the jet axis, suggesting that there is a deficit of X-ray 
emission near the radio lobes. The central gas density 
is ~ 4.6 x 10~ 3 cm~ 3 , temperature kT = 1.2 keV, and 
central pressure of 10 _12 Pa = 10~ n dyne cm~ 2 at a 
radius of about 0.5 arcminute. Consequently the pres- 
sure at a radius of few arcseconds from the nucleus could 
be a factor of 10 or so higher, or ~ 10~ 10 dyne cm~ 2 . 
A comparison between this number and the relation of 
Eqn. 1121 shows that the pressure in the X-ray emitting 
ISM is large enough to perturb the gas disk on reasonable 
timescales, so long as the isobars of X-ray emitting gas in 
the ISM are very elongated, structured, and anisotropi- 
c ally distributed. 

Qui Uen fc Bower! l)1999|) explored a simplistic model in 
which the ambient ISM has isobars elongated along the 
jet axis. In this case, the ISM exerts a torque on the disk, 
and gives rise to precession about the jet axis. However, 
if our model is correct, the disk in 3C 449 passes smoothly 
from the galaxy midplane at its outer edge to the plane 
perpendicular to the jet. If the jet causes a wind, the 
disk would be lifted out of the plane p erpendicular t o 
the jets, rather than pushed into it (e.g., lQuillenl l2001). 
Similarly, an azimuthal pressure variation across the disk 
may overcome the influence of the galaxy potential and 
perturb the disk at 100 pc, though not in ways capa- 
ble of creating the warp we predict e xists in 3C 449. 
As discussed bv lQuillen fc Bowerl l)1999)) . gas isobars ini- 
tially aligned with the jet should align with the galaxy 
equipotentials on a sound crossing time, of order the ro- 
tation period at 100 pc (a few Myr). Therefore the torque 
on the disk provided by static gas pressure in the ISM, 
no matter how strong or non-axisymmetric, could only 
cause the disk to precess and eventually settle into the 
galactic midplane, and not the plane perpendicular to 
the jets. The precessing disk may be damped via inter- 
nal kinematics (e.g., high viscosity, collisions, etc.), but 
the time associated with the damping would be of order 
the accretion timescale, which is too long to account for 
the warp. Even if the region is ram pressure dominated 
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and the gas in the disk is stirred near the sound speed, 
the azimuthal pressure gradient requ ired to warp the 
disk w ould be unrealistically large. The Quillen fe Bower 
( 1999) model is therefore incapable of explaining the ex- 
istence of the warp in 3C 449. That study, however, did 
not take into account the emergence of Chandra obser- 
vations revealing temperature fluctuations and buoyant 
bubbl es in the X-ray emi tting ISM surrounding an AGN 
(e.g., iFabian et all l2003|) . Such structure in the region 
is likely associated with the jet lobes, responsible for in- 
jecting large amounts of energy and momentum into the 
ambient medium. A drag force between jet-excited bub- 
bles and the di sk would not give rise to an instability 
(as explored bv iQuillenl 1200 lj) . and the warp would be 
damped. The inner disk could thus be pushed into the 
plane perpendicular to the jet. If the UV emission traces 
an excited cone of material above the disk and oriented 
along the jet axis, it may reveal the medium responsible 
for the torque that changed the alignment of the inner 
disk. 

5. SUMMARY AND DISCUSSION 

In this paper we have examined a near-infrared/ visible 
colormap of the radio galaxy 3C 449, finding a twist in 
the isochromes. This twist is unexpected, as a planar 
disk should have isochromes aligned with the disk major 
axis and exhibiting reflective symmetry about the mi- 
nor axis. The disk edge at a radius of ~ 600 pc has a 
major axis nearly (within 10°) coincident with the galac- 
tic isophotes (as seen at 1.6 /J.m), suggesting that the 
outer disk has settled to the midplane of the galaxy. The 
isophotes at 1.6 /im reveal no evidence for a recent merger 
or a strongly triaxial bulge. Consequently, the twist in 
the isochromes is most likely due to the disk itself. 

By integrating light through an absorbtive thin disk 
model, we find that we can account for the twist in the 
isochromes with a warped geometry. The model predicts 
rotation along the minor axis of the outer disk that is con- 
sistent with the sign of mean velocities measur e d from 
spectroscopic observations by iNoel-Storr et al.l {2003). 
In relation to the war p model, the pattern of UV 
and H a+PNIII emission (jAllen et aLll2002t iMartel et all 
2000), brighter to the north of the nucleus, could be emis- 
sion from above the disk, possibly in a wide cone oriented 
approximately along the jet axes. The UV emission in 
front of the disk and to the north of the nucleus would be 
easier to see than an opposing cone behind the disk and 
to the south. This suggests a possible explanation for the 
lack of UV emission south of the nucleus or associated 
with the disk on the western side. 

A warp in a disk with a sharp outer edge is not nec- 
essarily unexpected. For example, NGC 7626's con- 
tains a dust lane that ends abruptly, suggesting that 
it has an outer edge similar to that seen in 3C 449. 
Th e disk in NGC 7626 is cl early warped (see images 
bv lVerdoes Kleiin et al.lll999|) . and could resemble that 
of 3C 449 if viewed at a lower inclination. The ion- 
ized inner disk of NGC 6251 is tilted with respect to 
its outer dusty disk that also exhibits a sharp outer edge 
ijFerrarese fc Fordlll999|) . 

Our model for the warped disk of 3C 449 predicts that 
it is nearly parallel to the jet axis at a radius of 600 
pc, but twists and becomes nearly perpendicular to the 
jet at ~ 100 pc. We can find examples of other radio 



galaxies which might exhibit similar warps. For exam- 
ple, the disk in NGC 7626 is nearly edge-on, and twists 
such that the disk is nearly perp endicular to the jets 
near the nucleus (see images from iVerdoes Kleiin et alJ 
1999;. 3C 465, NGC 5127 and NGC 7052 contain disks 
with sharp outer edges. For these disks, as is true for 
that of 3C 449, the pattern absorption is almost trian- 
gular, suggesting that there is a twist in the isochromes. 
The jets in NGC 7052 are 45° from its disk, and the 
inner disk could be nearly perpendic ular to the jet (see 
images by Vcrdocs Kleiin et al. 1999). If the disks in 3C 
465 and NGC 5127 are twisted, their inner disks would 
be nearly perpendic ular to their jet axes (see images b y 
Ide Kofi et all l2000l and bv IVerdoes Kleiin et aljfl999|) . 
3C 236 contains an inner disk which is nearly perpendic- 
ular to its jets, whereas the o uter disk is tilted ab out 25° 
degrees from its outer disk l)de Koff et al.ll2000|) . M84 
3C 272.1) has a simila r inner and outer disk structure 
Quillen fc Bowerll9 99'l. and again the inner disk is more 
nearly perpendicular to the jet axis. NGC 6251 presents 
a counter example, as t he inner ionized disk is not ori- 
ented along the jet axis ijFerrarese fc Fordlll999|) . In the 
future, our procedure for modeling warped dusty disks 
could be used to study the disk geometries of these other 
objects. 

The notion that our warp model for 3C 449 accounts 
for jet/disk orthogonality is somewhat mysterious. Since 
the rotation period is shorter at smaller radii, an initially 
misaligned disk should first settle into the galaxy mid- 
plane at its innermost edge. However, the disk of 3C 449 
is settled into the symmetry plane at large radii but not 
in its inner regions. Such a morphology suggests that 
there is a force, aligned with or associated with the jets 
or active nucleus, that has pushed the disk preferentially 
at smaller radii. 

We have considered the possibility that a torque asso- 
ciated with a nuclear binary black hole could affect the 
disk orientation. We have also noted that radiation pres- 
sure from a central source could cause the disk to pre- 
cess, or that t he radiative warp instability proposed by 
IPringld l|1996j) is playing a role. We have shown that the 
timescales associated with these scenarios are too long 
to account for the warp in the disk of 3C 449. The pres- 
sure in the X-ray emitting ISM, however, is likely to be 
large enough that this medium could perturb the disk. 
A static model with elongated isobars would only cause 
the disk to precess, not be pushed into the plane orthog- 
onal to the jet axis. Such precession might be damped 
by processes internal to the disk (e.g., high viscosity, col- 
lisions, etc.), but the associated settling time would be 
of order the accretion timescale, which is again too long. 
Non-axisymmetric winds driven off the disk by evapora- 
tiv e, radiative, and MHD processes (like those discussed 
bv lMever fc Mever-Hofmeisterlll994T> are likely to play a 
role in the morphology of the disk, but not to a degree 
large enough to account for the warp. An interaction 
between the disk and the jet-excited ambient ISM (re- 
sulting in a drag force or inflow) could cause the disk 
to be pushed into the plane perpendicular to the jet, so 
long as there is sufficient turbulence and anisotropy in 
the X-ray emitting region. This suggests the intriguing 
possibility that studies of nuclear disk morphology can be 
used to probe past interactions between outflow-excited 
material and a comparitively cold gaseous disk. In the 
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case of 3C 449, the UV emission (visible only north of 
the nucleus) might help trace the intricacies of such an 
interaction. 

A close look at images of disks in radio galaxies re- 
veals not only warps but lopsidedness, sharp edges, and 
spiral features. Such observations are somewhat coun- 
terintuitive. Because of differential rotation, a lopsided 
disk should rapidly (on a few times the rotation period) 
become circular, a sharp edge should become smooth 
on a diffusion timescale (roughly the rotation period 
divided by the viscosity alpha parameter), and a non- 
self-gravitating disk should only exhibit transient spiral 
features. As such, the morphology of such disks sug- 
gests that they may have been recently perturbed, even 
truncated. Notably nonactive elliptical galaxies can con- 
tain disks which exhibit similar structures. It is possible 
that gaseous disks in elliptical galaxies are perturbed by 
present or past nuclear activity. Future high angular 
resolution studies could better probe this possible con- 
nection. 

Such a relationship between nuclear disk morphology 
and the ambient ISM can be considered a specialized 
manifestation of "feedback" from an AGN. Accreting 
black holes release enormous amounts of energy to their 
surroundings in various forms (through heating, radia- 
tion pressure, and induced kinetic motions, as discussed 



bv lBegelmar]|2004llF^bTan et alJl200^1 . It is not unrea- 
sonable, then, to interpret structurally perturbed nuclear 
disks as consequences of current and past interactions be- 
tween an AGN and its surrounding medium. 
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